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'nJ- ■ ABSTRACT 

\ We present a comprehensive study of 250,000 galaxies targeted by the Baryon Oscil- 

lation Spectroscopic Survey (BOSS) up to z « 0.7 with the specific goal of identifying 
CNJ ■ and characterising a population of galaxies that has evolved without significant merg- 

, ing. We compute a likelihood that each BOSS galaxy is a progenitor of the Luminous 

■ Red Galaxies (LRGs) sample, targeted by SDSS-I/II up z « 0.5, by using the fossil 

record of LRGs and their inferred star-formation histories, mctallicity histories and 
dust content. We determine merger rates, luminosity growth rates and the evolution 
. . of the large-scale clustering between the two surveys, and we investigate the effect of 

p% \ using different stellar population synthesis models in our conclusions. We demonstrate 

that our sample is slowly evolving (of the order of 2 ± 1.5% Gyr^^ by merging) by 
computing the change in weighted luminosity-per-galaxy between the two samples, 
and that this result is robust to our choice of stellar population models. Our conclu- 
sions refer to the bright and massive end of the galaxy population, with M^o.ss < —22, 
and > 10"-2Mq, corresponding roughly to 95% and 40% of the LRGs and BOSS 
galaxy populations, respectively. Our analysis further shows that any possible excess 
of flux in BOSS galaxies, when compared to LRGs, from potentially unresolved tar- 
gets at z w 0.55 must be less than 1% in the r"'^^— band (approximately equivalent 
to the band in the rest-frame of galaxies at z = 0.55). When weighting the BOSS 
galaxies based on the predicted properties of the LRGs, and restricting the analysis 
to the reddest BOSS galaxies, we find an evolution of the large-scale clustering that 
is consistent with dynamical passive evolution, assuming a standard cosmology. We 
conclude that our likelihoods give a weighted sample that is as clean and as close to 
passive evolution (in dynamical terms, i.e. no or negligible merging) as possible, and 
that is optimal for cosmological studies. 

Key words: galaxies: evolution - cosmology: observations - surveys 



© 0000 RAS 



2 Tojeiro et al. 



1 INTRODUCTION 

The Baryonic Oscillation Spectroscopic Survey (BOSS), 
part of the Sloan Digital Sky Survey III (SDSS-III), is an 
ambitious galaxy redshift survey which will determine the 
expansion rate of the Universe up to z « 0.7 by measur- 
ing the baryonic acoustic oscillations (BAO) and redshift- 
space distortions (RS D) in the galaxy power spectrum 
(|Eisenstein et al.]|201ll '). At the end of the five-year observ- 
ing program, BOSS will have mapped 1.5 million massive 
galaxies in 10,000 square degrees of sky, resulting in un- 
precedented volume and galaxy density. Forecasts indicate 
that BOSS will yield measurements of the redshift-distance 
relation dA{z) and of the Hubble parameter H{z) to 1% and 
1.8% at z — 0.35 and 1% and 1.7% at 2 = 0.55, respec tively 
(at the 1-sigma confidence level. .Eisenstein et al.ll20"lll ). Us- 
ing one third of the data, iReid et al.l (|2012h placed initial 
constraints on H{z = 0.57) and dA (Z = 0.57) at the 2.8% 
and 4.8% level respectively, whereas I Anderson et al.l (|2012l ) 
constrained Dv{z = 0.57) = [cz(l + zfd\H-^]^^^ to 1.7%. 
To achieve the survey's ambitious goals systematic uncer- 
tainties in the data, modelling and methodology must be 
kept to a minimum, and b e understood as best as possible 
(see IRoss et al.1 12OIII . |2012| for a study on data systematics 
in BOSS). 

A source of uncertainty in the modelling and measure- 
ment of the BAO is galaxy bias. Different populations of 
galaxies relate differently to the underlying matter density 
field, yielding different biases and often different scales that 
mark the regime over which a linear, deterministic and scale- 
invariant bias model is applicable. To a certain extent one 
can parametrise over this uncertainty, but nonetheless an 
interesting question remains concerning how much gain is 
possible if the bias modelling and evolution with redshift 
were well understood. 

The best candidate for a population of galaxies with 
a well understood bias evolution is a population that has 
been evolving with no or very little merg i ng: th e bias evo- 
lution is easily modelled using t he iFrvl (|l99m formalism 
(see also lTegmark fc Peebleil998l '). Massive red galaxies are 
the prime candidates for such a population - they are com- 
posed mostly of old stellar populations (e.g. iMaraston et al.l 
I2OO9I ). and their growth via merging since a redshift of two 
has been con strained to be small (< 10%) even if strictly 
non-zero (e.g. I Wake et a l. 2008). Halo-modelling analyses of 
massive red galaxies have repeatedly revealed a highly bi- 
ased population (6 « 2) with a low satellite fraction (5 t o 
10% of galaxies ar e satellites) - see e.g. Zehavi et al.l (l2"005al) : 
IWake et al.1 (|2008l ): IZheng et al.1 (|2009[ ): IWhite et all \2Qm 
- confirming their suitability for cosmological studies. De- 
parture from a pure passive dynamical evolution history 
has bee n shown to have a dependence on luminosit y and 
colour (|Toieiro fc Percivall [2OI0I : iToieiro erall l201ll '). and 
this opens up the possibility of weighting galaxies appro- 
priately, so as to maximise the contribution of those that 
are more likely to have been passively evolving, and min- 
imise the contribution of those that are less likely to have 
done so. The SDSS-I/II survey targeted LRGs using a mix 
of colour and luminosity selection cuts such as to follow the 
evolution of a passively evolving population of stars. BOSS 
targeting, however, is much less restrictive in terms of lumi- 
nosity and colour (especially a,t z > 0.45 - see Section (2]). 



It is therefore not true that one population is automatically 
composed of the evolved products from the other. 

One of the goals of this paper is to identify, in BOSS, the 
most likely progenitors of lower redshift SDSS-I/II LRGs, 
and design a set of weights that allow a selection of the 
galaxies that are linked by the same evolutionary history. 

Our other major goal is to place quantitative constraints 
on the formation and recent evolution of present-day lumi- 
nous red galaxies, which in broad terms constitute a sub- 
set (at large luminosities, or stellar masses) of what are 
typically called early-type galaxies (ETGs). Efforts towards 
understanding ETGs and their evolution can be split into 
two categories: those that focus on their stellar content, 
and those that primarily aim to constrain their dynamical 
evolution, or merging history. Studies have been performed 
based on (see also references within): the mass or luminosity 
funct i on of central galax ie s JWake et al.l 20061 : iBrown et al.l 
I2OO7I : iFaber et al] 120071: [Cool et al.l 12008^. and of their 
\20m~ 



satellites (iTal et al 



col our- m agnitude diagram 



JCool et al.1 12OO6I: iBernardi et al.] I2OIII'): photora etry SED 
fitting ( Kavirai et al. 20091 : Maraston et alj|2009|): absorp- 
tion li ne fitting to individua l galax i es' spectra jTrager et al.l 
I2OOOI : iThomas et al] |2005|. I2OI0I: ICarson fc Nichoj I2OI0I) 



or to stacked spectra (jEisenstein et al.l 1200 



Graves et al 



200d: IZhu et al.ll20ld '): full spectral fi t ting jjimenez et al 
2007); close-pair counts (Bell et al' '2006^: 'Bundy et al 
2009) and clustering (Zehavi et al. 2005a; Shcth ct al. 20oi; 



Masiedi et al.l I2OO6I : IConrov et aLjl2007l: IWhite et al 
Brown et alTboOSi : iMasiedi et al.l boosi: I Wake et al. 



2007; 



20081 : 



Toieiro fc Percivall [20T0I : iDe Propris et al.1 l2010h . There is 
general agreement in the overall picture: ETGs constitute a 
uniform population of galaxies; are dominated by old and 
metal rich stellar populations; their mean ages (either mass- 
or light-weighted) decrease with luminosity; and the most 
luminous occupy the more dense environments. There is, 
however, an increasing amount of evidence pointing towards 
some amount o f recent star formation i n intermediate-mass 
ETGs (see e.g. ISchawinski et al.|[2007l : iKavirai et al.ll2007l : 



ij^iLis (see e.g. iscnawmski et ai.M:^UUYl : iKavirai et ai.M:^UUYl 
ISalim fc Richll2010l ). This amount of star formation is not ir 
con flict with the h i erarch ical model of structure formation, 
and iKavirai et al.l l|2010 ''). through evidence coming from 
small morphological disruptions in early-type galaxies, ar- 
gue that it can be explained from the contributions from 
minor-mergers. 

On the clustering side, halo modelling is rapidly be- 



formation (see e.g. Zehavi et al.ll2005b|; Zheng et al. 


2007 


Skibba fc Sheth 2009|; Skibba 2009: Ross fc Brunnei 


200S 


Zheng et all 120091: Ross et al.1 12010: iTinker & Wetzel 


201c 


Wake et alj I2OIII. and references within). It is a power- 



ful approach that connects galaxies with the dark matter 
halos in which they reside, and which describes the dis- 
tribution of a population of galaxies in terms of centrals 
and satellites, as well as their relative ratio, as a func- 
tion of hal o mass (which is we ll correlate d with luminos- 
ity, see e.g. lSwanson et all 20081; [C rcsswcU fc Percivalll2009l : 
IRoss et alj|201ll '). E.g., IZheng et al.. (,20070 use luminosity 
dependent galaxy clustering at different epochs and the ex- 
pected growth of dark matter halos to infer a growth due 
to star formation between 2 = 1 and the present day, after 
roughly taking into account growth due to the merging of 
centrals and satellites. This type of description of galaxy as- 
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sembly can then be directly co mpared to predicti ons from 
semi-analytical simulations (see lZehavi et al1l2012l ). 

More specifically, the dynamical passive model can be 
directly tested by a halo model type of analyses. By per- 
forming HOD modelling at two different redshifts, one can 
evolve the best-fit halo model fitted at one redshift to an- 
other, assuming passive evolution. Comparison of the best- 
fit halo models provides insight about the dynamical evolu- 
tion of the sample, particularly in terms of satellite accre- 
tion and disruption. For most of the samples chosen, anal- 
yses show that a purely passive model would predict too 



must merge or be disrupted (see e.g. 
White ct al 2007; Zheng ct al. 2007; 



'Conrov et al.l 


2007 


Brown et al.l 


2008 



Wake et al]|2008l : ISeo et ai.ll2008l '). Measurements of merger 
ra tes of massive galaxies vary significantly (see Table 4 
in lToieiro fc Perciva]|[2O10l for a summary), but luminosity 
growth via merging seems confined to something between 
3-20% since z x 1. 

It seems increasingly likely that the assembly history 
of massive galaxies is inexorably linked to the existence of 
intra-cluster light (ICL) - a diffuse and scattered stellar com- 
ponent that can account for 10-50% of the stellar mass in 
clusters (see e.g. iFeldmeier et al.l |2004| : iMihos et all l2005l : 
iPurcell et aLll2007l : lYang et all 120091 '). A likely mechanism 
of its formation is the disruption of satell ite galaxies when 
halos merge (see e.g.lCo nroy ct al. 2007; P urcell et al]|2007l : 
I White et al.ll2007l : lYang et al.ll2009l and discussions therein). 
A lack of conservation of light, or stellar mass, in galaxy 
mergers has implications for the interpretation of the evo- 
lution of the luminosity function and inferred merger histo- 
ries. The fraction of light lost by a merging satellite to the 
intra-cluster medium remains largely unconstrained, with 
estimates at the large halo mass end from the studies cited 
above varying between 15% and 80%. In the present work 
we make no explicit allowances for the loss of light to the 
ICM when two galaxies merge, but we will argue that our 
results are robust to this eff'ect, within the limitations of the 
models and data. 

In the work presented here we approach the problem 
of galaxy assembly from a ne w direction, opposite in ethos 
to that of IZheng et al.] (|2007l ). We will use state of the art 
modelling of the stellar evolution of a sample of galaxies 
to directly quantify growth from star formation, and from 
that infer a galaxy-merger history. We compute a model for 
the stellar evolution of SDSS-II LRGs by decomposing their 
spectra into a series of star-formation and metallicity histo- 
ries, as well as dust content. This allows us to make predic- 
tions of their colour and magnitudes at any redshift. This 
information, when combined with the target selection infor- 
mation for BOSS galaxies, constrains the regions in colour 
and magnitude space in BOSS within which progenitors of 
LRGs are more likely to reside. We then compute a set 
of weights that depend on the predicted evolution of each 
galaxy across the two surveys, and up-weight the objects 
that are more likely to be in both samples. The analysis 
we present depends on underlying assumptions about stel- 
lar evolution, initial mass functions and dust modelling. We 
perform the full analysis using two different sets of assump- 
tions, so as to give the reader an idea of the dependence our 
final results on this type of uncertainty. 

Isolating the likely progenitors of LRGs in BOSS is in 



itself no test of the mer ging history of the sample . Follo wing 
on from our analyses in lToieiro &: Percivall ([2013, [20l3), we 
test the evolution in the number and luminosity density of 
the galaxies between LRGs and BOSS, as a way to mea- 
sure the amount of merging or luminosity growth between 
the two redshift surveys. We also use a luminosity- weighted 
two-point correlation function to further test the dynamical 
passive hypothesis - weighting the galaxies by luminosity 
produces a clustering statistic that, on large-scales, is less 
sensitive to galaxies within the sample merging. 

This paper is organised as follows: in Section [2] we de- 
scribe our two data sets, including targeting; in Section [3] 
we explain how we compute a stellar evolution model that 
describes the stellar evolution of all galaxies and spans a 
redshift range between 0.23 and 0.7; in Section U we use 
this stellar evolution model to compute a set of weights that 
allows us to construct optimal samples of galaxies at dif- 
ferent redshifts and explore the evolution of LRGs in the 
BOSS volume; in Section [5] we compute merger-rates and 
average luminosity growth across the samples and, in Sec- 
tion |S1 we compute the large-scale clustering of each of our 
samples and compare to predictions from a purely passive 
model. Finally we discuss and summarise our conclusions 
in Section [7] Where required we assume a fiat A cold dark 
matter (LCDM) cosmology with Qm = 0.266, = 0.734 
and Ho = 71 kms~^Mpc"^ 



2 DATA 

The Sloan Digital Sky Survey (SDSS) has imaged over 
one quarter of the sky using a dedicated 2.5m telescope in 
Apache Point, New Mexico l|Gunn et al.ll2006l) . For details 
on the hard ware, software an d data- reduction see lYork et al.l 
(|2000l ) and IStoughton et all l|2002l ). In summ ary, the sur- 
vej^was carried out on a mosaic CCD camera l|Gunn et al.l 
[19981 and an auxiliary 0.5m telescope for photometric cali- 
bration. Photometry was taken in five bands: u, g, r, i and z 
l|Fukugita et al.ll 19961 ). and magnitudes corrected for Galac- 
tic extinction using the dust maps of Schlegel et al.l I 



BOSS, a part of the SDSS-III survey (|Eisenstein et al 



199S 



2011 



has mapped an additional 5, 200 square degrees of southern 
galactic sky, increasing the total imaging SDSS footprint to 
nearly 14, 500 square degrees, or just over one third of the 
celestial sphere. All of the imaging was re-processed and re - 
leased as part of SDSS Data Release 8 l|Aihara et al.ll201ll ). 

In SDSS-I/II, Luminous Red Galaxies (LRGs) were se- 
lected for spectroscopi c follow-up acco r ding to the target 
algorithm described in lEisenstein et aU l|200l[ ), designed to 
follow a passive stellar population in colour and apparent 
magnitude space. In this paper we analyse the latest SDSS 
LRG spectroscopic sample (Data Release 7. lAbazaiian et al] 
l2009t ). which includes around 180,000 objects with a spec- 
troscopic footprint of nearly 8000 sq. degrees and a redshift 
range 0.15 < z < 0.5. In SDSS-III, the BOSS target selec- 
tion extends the SDSS-I/II algorithm to target fainter and 
bluer galaxies in order to achieve a galaxy number density 
of 3 X 10"* h^ Mpc"^ and increase the redshift range out to 
z « 0.7. The spectroscopic footprint of the BOSS data used 
in this sample covers almost 3500 sq. degrees of sky, and cor- 
responds to the upcoming Data Release 9, which will mark 
the first spectroscopic data release of BOSS. 
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The targeting algorithms make use of five different def- 
initions of magnitudes as follows: 

• SPSS uber-calibra ted model magnitudes 

l|Padmanabhan et al.l |2008| ). computed using either an 
exponential or a DeVaucouleurs light profile fit to the 
r-band only, denoted here with the rnod subscript; 

• cmodel magnitudes, computed using the best-fit linear 
combination of an exponential with a DeVaucouleurs light 
profile fit to each photometric band independently, and de- 
noted here with the subscript cmod', 

• point-spread function (PSF) magnitudes, denoted with 
a psf subscript, and computed by fitting a PSF model to the 
galaxy; 

• petrosian magnitudes, computed from the petrosian 
fiux (the fiux measured within twice the Petrosian radius, in 
turn defined using the surface brig htness of the galaxy, see 
|Petrosiajilll976l : IStrauss et al.ll2002l ). and denoted here by a 
subscript p; and finally 

• fibre magnitudes, computed within a 2 arcsec aperture, 
and denoted by a fib2 subscript. 

In SDSS-I/II, redshift LRGs were selected using two dif- 
ferent algorithms. Cut-I predominantly but not exclusively 
targeted lower redshift galaxies {z < 0.43) using the follow- 
ing selection criteria: 



rp < 13.H-C|| (1) 

T-p < 19.2 (2) 

c_L < 0.2 (3) 

fir,p < 24.2 mag arcsec'^ (4) 

rpsf - rmodei > 0.3, (5) 
where the two colours, cy and c± are defined as 

C|| = 0.7(3 -r) + 1.2[(r-i)- 0.18] (6) 

c± = (r- - i) - (5 - r)/4 - 0.18. (7) 



Model magnitudes are used for the colour cuts, and pet- 
rosian magnitudes for the apparent magnitude and surface 
brightness constraints. Note that whereas petrosian mag- 
nitudes naturally fail to account for fiux outside twice the 
petrosian radii, and whe reas this fraction var ies as a function 
of galaxy type (see e.g. iGraham et al.l [20051 ). here they are 
simply used to define a sample of galaxies. When comput- 
ing luminosity densities we always use cmodel magnitudes. 
Cut II mostly but not exclusively targets LRGs at 2 > 0.4 
following: 



T-p < 19.5 (8) 

c_L > 0.45 - (.g - r)/6 (9) 

(g-r) > 1.3 -h 0.35(r - i) (10) 

^r,p < 24.2 mag arcsec'^ (11) 

Tpsf - Tmodel > 0.5. (12) 



Two separate algorithms are necessary as the passive 
stellar population turns sharply inai; — rvsr — i colour 
plane, when the 4000A break moves through the filters. 

In SDSS-III, galaxies ai z < 0.43 are predominantly 
but not exclusively targeted by the LOZ selection algorithm. 



akin to Cut I above, but extended to fainter magnitudes. A 
LOZ galaxy must pass the following: 

r-cmod < 13.6 + C||/0.3, (13) 
|cx| < 0.2, (14) 
16 < Tcmod < 19.6, (15) 

where the two auxiliary colours cy and c± are defined as for 
SDSS-I/II above. 

Galaxies at 2 > 0.43 are predominantly but not exclu- 
sively targeted by the CMASS selection algorithm, which 
extends the Cut II above by targeting both fainter and bluer 
galaxies. A CMASS galaxy must pass the following criteria: 



17.5 < iemod < 19.9, (16) 

< 2, (17) 
d_L > 0.55, (18) 

< 21.5, (19) 
icmod < 19.86 -fl.6(dj,erp- 0.8), (20) 

where the auxiliary colour d± is defined as 

d± = rmod — imod — {gmod ~ rmod)/8.0. (21) 

CMASS objects must also pass the following star-galaxy 
separation cuts: 

ipsf-hnod > 0.2 + 0.2(20.0- w), (22) 

Zpsf - Z,nod > 9.125 - OAGZmod, (23) 

unless they also pass the LOZ cuts. 



The CMASS selection algorithm was designed to loosely 
follow a constant stellar mass limit and, unlike Cut-II in 
SDSS-II, it does not exclusively target red objects. There- 
fore, whereas both the LRG and CMASS samples are colour- 
selected, CMASS is a significantly more complete sample 
than the LRGs, especially at the bright end. In this pa- 
per we will split our data into two distinct redshift slices, 
with our lower redshifts slice ranging between 0.23 < 2 < 
0.45 and our higher redshift slice between 0.45 < 2 < 
0.7. Our low redshift slice consists exclusively of SDSS- 
I/II LRGs (Cut-I and Cut-II) and contains approximately 
89,000 galaxies, and our high redshift slice consists exclu- 
sively of SDSS-III CMASS galaxies, with over 250,000 ob- 
jects. The low redshift cut-off is motivated by our previous 
analysis of the LRGs that indicates t he sample is signifi- 
cantly contaminated at lower redshifts (|Toieiro et al.|[201ll : 
iToieiro fc Percivalllioill ') . We do not make use of LOZ galax- 
ies for the main analysis presented in this paper, mainly due 
to the fact that the volume and number density sampled by 
LOZ currently lags behind that of the CMASS due to prob- 
lems in target selection at the beginning of the observing 
run. We use LOZ galaxies only in Section [5.41 when investi- 
gating potentially unresolved targets in CMASS. The n{z) 
distribution of our two samples is shown in Fig. [T] 



3 THE STELLAR POPULATION MODELLING 

We use the 124 s tellar evolution models computed in 
iToieiro et al.l (|201ll ) by stacking LRG spectra according to 
their luminosity, redshift and colour, and subseque ntly anal- 
ysed them with VESPA (jToieiro et al.ll2007l . [2OO9I ) to obtain 
detailed star-formation histories as a function of lookback 
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Figure 1. Number density as a function of redshift for the LRG 
(red) and the CMASS (black) samples. The dashed line a.t z = 
0.45 shows our chosen hard boundary between the two surveys - 
we do not use any LRGs with z > 0.45 nor any CMASS galaxies 
with z < 0.45 



time. VESPA fits a linear combination of stellar populations 
of different ages and metallicities, modulated by a dust ex- 
tinction, to the stacked optical spectra. Each star-formation 
history can then be translated into a detailed evolution of 
any magnitude and colour with cosmic time. We have made 
no changes to these publicly available models other than 
increasing the sampling in redshift, to provide better re- 
solved colour and magnitude evolutior|j. We consider the 
solutions obtained with two sets of stellar population mod- 
el s: the Flexible Stellar Pop ulation Synthesis (FSPS ) models 
of lConrov et~aLl ([2009.) and lConrov fc GunnI li2010h. and the 
stellar population models of iMaraston &: Strombaci3 (|201ll) 
(Mil) - we refer the reader to Section 4 of Toieiro et al.l 
l|) for detailed information on the differences and sim- 



(|201 

ilarities between the two sets of assumptions, and we note 
that the most significant difference arises fro m the stellar 
evolut ion tracko One of the main results in iToieiro et al.l 
(|201l[) is that, even though FSPS and Mil provide star- 
formation histories that have very similar mass-weighted 
ages that decrease with luminosity, in the Mil case this 
is due to the presence of a population of stars of young to 
intermediate ages (1-3 Gyr), whilst in the FSPS case this is 



^ The models from Tojeiro et al. (2011) are available at 
|http: //www. leg .port . ac.uk/- tojeiror/lrg.evolution/ " 

^ Briefly, notable differences lie in the choice of the shape 
of the initial mass function (IMF), isochrone tracks and 
stellar lib raries. I n the case of Mil, we use a combina- 
tion of a iKroupal ('2001) IMF, the MILES stell ar library of 
ISanchez-Blazquez et al. (20 06) and isochrones from lCassisi et al] 
and ISchaller et ah I ll 19921') combined with the fuel con- 
sumption approach of iRenzini fc Buzzonil ( Il986t) for post Main- 
Se quence p h ases. For FSPS models we use a combination of 
a IChabrieJ ||2003|') IM F with a MILES stella r library and 
the P adova isochrones of lMarigo fc Girardil ( |2007| ) : [Marigo et al.l 
1I2OO8I ') 



due to a slightly younger main burst of star formation, which 
extends to lower ages with decreasing luminosity. These dif- 
ferences in the star-formation histories will have an impact 
on the results, and we will compare results obtained using 
both models throughout the paper. In Section [3. 1.1 1 we de- 
scribe the star-formation histories recovered with both sets 
of models in detail. 

In Fig. [2] we show the g — r and r — i colours predicted 
by the fits to LRGs based on the different models (red dots 
for FSPS and blue for Mil) and how they compare to the 
colours of observed CMASS galaxies in four redshift ranges 
(blue contours). The locus of the models traces the locus 
of the observed galaxies remarkably well. Furthermore, the 
FSPS models predict a tendency to have bluer colours with 
increased redshift, and that is tentatively matched by the 
data. Mil models follow broadly the same trend, with the 
main differences seen at z = 0.55, where Mil models predict 
significantly bluer galaxies (some models predict a crossing 
of the d± cut). 



3.1 The composite model 

The clear advantage of our set of models is that it gives 
a data-driven grasp on the stochasticity of the population 
properties. We do not need to assume all targeted LRGs 
are the same and natural scatter in the colours - given by 
changes in metallicity and star-formation rate - can be triv- 
ially accounted for. We are limited in the sense that we 
can only predict the evolution of any galaxy to redshifts 
greater than the one it is observed at; this is because the 
fossil record can only hold information on the past history 
of a galaxy. Evolving a galaxy forward requires assumptions 
about any subsequent star- formation, or lack of it. In order 
to match the samples we need a stellar evolution model that 
spans the redshift range of both samples combined, and that 
we can use to evolve any galaxy to any redshift with min- 
imal assumptions about their stellar evolution. We choose 
an approach where we compute a single weighted composite 
model that spans the redshift range of the sample, based 
on the 124 individual stacks. At each redshift we compute 
a mean spectrum, weighted by the number of galaxies that 
make a prediction for that particular redshift (i.e. observed 
at z < z'). From this mean spectrum we compute a new 
set of magnitude and colours, that define what we will call 
our composite model. The k-l-E corrections of the compos- 
ite model are the weighted means of the individual k-l-E 
corrections - this composite model is therefore our best esti- 
mate of the overall average colour and magnitude evolution 
of the full LRG sample. Note that this approach is formally 
the equivalent to taking the weighted mean of the 124 star 
formation histories for each stack, and using that weighted 
star formation history to recover the composite spectrum 
and correspondent models. 

We show the colour evolution of our model in Fig.|3l and 
the K-|-e corrections in Fig.|?](red for FSPS models, and blue 
for Mil). These models are used to describe all galaxies in 
the study: LRGs and CMASS galaxies alike. For complete- 
ness in Section [5.5l we briefiy discuss the im pact of using the 
strictl y passive stellar evolution model of IMaraston et al.l 
l|2009l ), or the full range of 124 individual models, on our 
results. 
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Figure 2. The observed colour evolution of CMASS galaxies contrasted with the predicted colour evolution of LRGs at CMASS redshifts. 
In each panel the black contours show the number density of the full LRG sample in the g — r vs r — i plane. The blue contours show 
the number density of CMASS galaxies for a given redshift range (given for each panel). The red dots show the predicted colours of the 
LRGs at the same redshifts given by the FSPS models, and the blue dots show the predicted colours using the Mil models. The different 
dots correspond to the prediction of LRGs of different luminosity, colour and redshift. The solid red line shows the dj_ = 0.55 cut for 
reference. 



3.1.1 Physical model 

As mentioned in the previous section, VESPA solutions with 
the two different stellar population models give physical 
models for the galaxies that are qualitatively different, espe- 
cially for LRGs at 2: < 0.25. FSPS produces a model that is 
nearly completely passive, with less than 2 — 3% by mass in 
stars that are younger than 3 Gyrs. Mil gives a model that 
sees over 90% of the stellar mass formed over 12 Gyrs ago 
(for a galaxy at z — 0), but which often puts a non- negligible 
amount of stars at ages of 1-3 Gyrs (up to 10% in mass). 
This generates more scatter in the blue points in Fig. [2] and, 
as a direct consequence, a larger scatter in Figs. O and 2] 

Small but non-negligible amounts of star-formation act 
to steepen the luminosity evolution (given by the k-fE cor- 
rections), as the galaxy effectively 'loses' stars as we step 
back in redshift. More generally, a change in k-l-E corrections 
can also arise from different assumptions in the stellar evo- 
lution models, or from a different slope - or an evolving slope 
- of th e initial mass function (IMF). In iToieiro fc Percivall 
l|20ldl we investigated the effects of an added redshift- 
dependent term to the k+E corrections, being motivated at 
the time by uncertainties in the slope of the IMF. Here we 
will perform no such investigation, but having two models 



with two different slopes for the k-f E corrections provides 
an estimate of the impact of this uncertainty on our final 
results. 

Both stellar population models give a constant metallic- 
ity with redshift, although Mil solutions are slightly more 
metal rich at Z ^ 0.03, whereas FSPS prefers a solution 
with Z ^ 0.025. 

Finally, the dust content is very similar in both cases 
- extinction increases with decreasing luminosity, increas- 
ing redshift and increasing r — i colour, varying between 
TV = 0.2 and tv = 0.8. Here tv is the optical depth at 
A — 550 0A and the du s t ext inction is modelled accord- 
ing to a Chariot fc Falil (|200(]| ) mixed-slab geometry (see 
iToieiro et al.ll201ll for lull details). The weighted average, 
and the effective extinction for the composite model, is 
TV ~ 0.4 - 0.5. 



3.2 K+e corrections 

We follow closely the procedure of IToieiro fc Percival 
(2010. ). which we summarise here for completeness. 

Our composite model provides L\{tage), the luminosity 
per unit wavelength of a galaxy of age tage- We K-fe-correct 
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Figure 3. The composite stellar evolution model, computed ac- 
cording to the procedure in Section [XT] In all panels the shaded 
contours show the number density of LRGs (at z < 0.45 and 
on the bottom half of the last plot) and CMASS galaxies (at 
2 > 0.45 and on the top half of the last plot). The red (blue) 
solid line shows our composite stellar model obtained using the 
FSPS (Mil) VESPA star-formation histories. It is a weighted av- 
erage of the models shown in Fig. [2] The error bars show the Icr 
dispersion of the models shown in Fig. [2] in each redshift bin. For 
reference, the yellow l ine shows the LRG purely passive model of 
iMaraston et al ] l|2009l ') - see Section [SXTI 



all galaxies to a common redshift of Zc — 0.55, and calculate 
corrected absolute magnitudes in filters shifted to Zc — 0.55 
as 



MrO.55 = rcmod 

with 

Ke{z, Zc) = 
= -2.5 log 



Slogi 



lOpc 



Ke{z,zc), (24) 
(25) 



0.4 0.5 
redshift 

Figure 4. K-fe corrections in the ro.ss-band (triangles) and in 
the io.55-band (asterisks). The red lines refer to the FSPS models, 
and the blue lines to the Mil models. The error bars show the Icr 
scatter around the mean from the 124 individual stacks. These 
corrections allow us to compute the evolved absolute magnitude 
0.7 of any galaxy at 2 = 0.55, in the two shifted filters (therefore 
for galaxies at 2 = 0.55 this correction is fixed and independent 
of their spectra or modelling). The corrections in the ro.55 band 
are steeper because it traces the 4000A break at these redshifts - 
see Fig. [5] The scatter in the Mil k-|-E corrections is larger, as 
these models predict stochastic events of star-formation at young 
to intermediate ages in some of the stacks. For re ference, the yel- 
low line shows the LRG purely passive model of lMaraston et al] 
(2009), as a dashed line for the jQ.ss-band and as a solid line for 
the ro.55-band - see Section l5.5.1l 



Ao is in the observed frame and Ae in the emitted frame. 
T\ is the SDSS's r-band filter response, and L\{z) the lu- 
minosity density of a galaxy at redshift z, given the fiducial 
model. We also compute Mio.55, using exactly the same pro- 
cedure on the i-band. Note that for a galaxy at z = 0.55, 
the K-|-e correction is independent of the observed or mod- 
elled spectrum and equals — 2.51og]^Q (^^Tz^- By choosing 
Zc = 0.55, roughly the peak of the redshift distribution of 
CMASS galaxies, we minimise the effect of the modelling 
on CMASS galaxies. The other option would have been to 
k-fe correct to median redshift of LRGs. However, as the 
composite stellar population model is based on the spec- 
tra of LRGs, its predictions must be at least as robust for 
LRGs as for CMASS galaxies, if not more so. Therefore, our 
procedure is the more robust approach. We show the k-l-e 
correction in the r— and i— bands in Fig. [4] For reference, 
in Fig. [S] we show the expected observed-frame spectrum of 
a typical galaxy in the sample at Zc = 0.55, along side the 
three broadband filters used in this paper. 



3.3 Comparing CMASS galaxies and LRGs 

We can use the K-|-e-corrected absolute magnitudes to 
broadly characterise the two samples. Fig. [S] shows a simple 
comparison of the magnitude distributions for both sam- 
ples and their evolution with redshift computed for both 



© 0000 RAS, MNRAS 000, 000-000 



8 Tojeiro et al. 




3000 4000 5000 6000 7000 8000 9000 10000 
wavelength [AA] 

Figure 5. The expected observed spectrum of a typical galaxy 
in the sample at z = 0.55 (black). The three broadband filters 
used for target selection are overplotted: g— band in blue, r— band 
in green and i— band in red. For reference, we show in grey the 
expected observed spectrum of a galaxy at ^ = 0.3. 



the r— and the i— band. Once again we show the results 
for the FSPS model in red and for the Mil model in blue. 
Here the only model differences come through the K+e cor- 
rections, with the different slopes between models (shown 
in Fig. [4]) naturally giving different k+E corrected abso- 
lute magnitudes. Mil shows a steeper slope with respect 
to FSPS, with the crossing point at Zc = 0.55. So for a 
galaxy a.t z < 0.55 Mil will predict a fainter k+E cor- 
rected magnitude at z = 0.55. Conversely, for a galaxy at 
z > 0.55, Mil will predict a brighter k+E corrected mag- 
nitude at 2 = 0.55. By construction, the magnitudes for 
galaxies sitting a,t z — 0.55 will match for both models due 
to our choice of filters. The top panel of Fig. |6] shows the 
effect of having different slopes for the k-l-E corrections - 
for LRGs this is about 0.3 magnitudes in the ro.ss-band; for 
CMASS galaxies it is much smaller, at less than 0.1 mag- 
nitudes. These values are roughly halved for the io.ss-band. 
The bottom two panels of Fig. |S] show the evolution of the 
corrected magnitudes with redshift (solid contours for FSPS 
and line contours for Mil). As expected, we see a steeper 
evolution with redshift using the Mil contours. 

Fig. [7] displays colour-magnitude relations. Here we 
show only the results using FSPS models as the results are 
similar in both cases. The CMASS sample has a broader 
range in absolute magnitude and colour than the LRG sam- 
ple, as expected given the larger number density. The clear 
trend seen between rest-frame colour and Mro.55 is explained 
simply by target selection. To help make this point we show 
the expected evolution of the colour-magnitude relation of 
an object at the faint end of the survey (cmodel — 19.9 at 
z = 0.45) and an observed colour of r — i — 0.8, between 
z = 0.25 and z = 0.7 - this is the red line in both plots. Any 
object to the faint side of the red lines would fail the mag- 
nitude cuts in the i— band of the CMASS algorithm. This 
gives an obvious artefact when plotting Mro.55 vs colour, 
where upon the CMASS selection does not select faint blue 



galaxies. The bright end slope is a consequence of volume 
effects, coupled with the slope of a typical galaxy spectra. 



4 SAMPLE MATCHING 

We now construct galaxy samples at high and low redshift 
that are coeval according to our composite stellar evolution 
models. We continue to closely follow the methodology of 
iToieiro fc Percivall l|2010D . which we summarise below. We 
have to take into account three redshift-dependent effects: 

(i) the intrinsic evolution of the colour and brightness of 
the galaxies; 

(ii) the varying errors on galaxy colour measurements; 
and 

(iii) the varying survey selection function. 

Our correction for (i) is given by our composite stel- 
lar evolution model. We include an evolvi n g colo ur scatter 
term to allow for (ii). iToieiro fc Percivall ( 2010l ) used the 
populati on scatter around the ste llar evolution model with 
redshift. Plbieiro fc Percivall |2oi3) updated this term to be 
based on the evolution of photometric errors as a function of 
apparent magnitudes, which were modelled as a function of 
redshift - see their Section 3. The motivation was two fold: 
firstly the photometric errors are driven principally by the 
apparent magnitude of an object, rather than its redshift; 
and secondly this is less dependent on choice of stellar evo- 
lution modelling. We adopt this approach here. For (iii) we 
construct a set of weights that assures a given population 
of galaxies - in terms of colour and absolute magnitude - is 
given the same weight in the high and low redshift samples, 
as described in the next section. 



4.1 Weighting scheme 

We use the weighting scheme of IToieiro fc Percivall l|2010l ). 
which keeps the total weight of each galaxy population the 
same in different redshift slices. 

Suppose an LRG, gA, is faint and therefore can only be 
seen in a small fraction of the CMASS volume, fv, but can 
be seen in the full LRG volume. Then our weighting scheme 
will give gA a weight that is equal to fv- Consider now a 
faint CMASS galaxy, that is observed in fv, and whose mag- 
nitude and colour evolution matches those predicted for gA- 
This galaxy will by definition also only be observed in a frac- 
tion fv of the CMASS volume. Our weighting scheme gives 
(?s a weight on unity. Note this is the opposite approach to 
the traditional Vmax weight, which would up-weight gs by 
1/ fv and give gA a weight of unity. 

Explicitly, for an LRG in a volume Vlrg we calculate 



Vlrg 



LRG 



LRG 



CMASS 



Vlrg Vcmass 



(26) 



and similarly for a CMASS galaxy, in a volume Vcmass'- 



Vr 



match, i — 



Vcmass 

yCMASS 
max.i 



* m 



• m 



Vlrg Vca 



(27) 



where Knax.i is the volume a galaxy i would have been ob- 
served in either survey, according to the full target selection 
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Figure 6. Comparing K+e corrected magnitudes in SDSS-I/II LRGs and BOSS CMASS galaxies. Top: the distribution of absolute 
magnitudes for LRGs (dashed lines) and CMASS galaxies (solid lines). The different colours show the results from using different stellar 
population models, with FSPS in red and Mil in blue. The two panels show the magnitude computed either in the rest-frame r— or 
i— band. Bottom: the absolute-magnitude with redshift on both samples. FSPS results are shown in the solid contours, and Mil in the 
line contours. The samples are split at ^ = 0.45; we do not use any LRGs with z > 0.45 nor any CMASS galaxies with z < 0.45. These 
plots show clearly the reach to fainter magnitudes of the CMASS sample. See main text for a discussion on the effect of the stellar 
population models. 



cuts and the evolution of its colour and magnitude, as given 
by the composite model. 

Where the traditional Knax estimator would up-weight 
galaxies only visible in a fraction of the volume they were 
observed in, we instead give these galaxies a weight of unity 
and down-weight the corresponding galaxies with the same 
properties observed in the other volume. 

The interpretation of the Knatch weight is different than 
that of the traditional l^max weighting. Whereas the lat- 
ter gives us the means to correct for incompleteness and 
yields true space densities, the former must be interpreted 
as a weighting scheme rather than a completeness correc- 
tion. I.e., Vmatch weighted number and luminosity densities 
are still potentially volume incomplete, but the populations 
are weighted in such a way that they are equally represented 
at both redshifts. We can compare the distribution of total 
weighted luminosity for the two slices, but we cannot inter- 
pret these functions as giving the true luminosity density. 

The advantage of this weighting scheme is that we sam- 
ple different populations equally based on volume, and there- 
fore obtain a weighted population such that galaxies ob- 
served throughout a large volume are up-weighted. It also 
implicitly checks that we are only using populations that 



exist in bot h samples, withou t having to do such a test ex- 
plicitly (e.g. lWake et al.ll2006l ). 



4.2 The progenitors of LRGs 

A large value of Vmatch (Vmatch varics between and 1) in- 
dicates that a galaxy belongs to a population that can be 
observed across a large fraction of both surveys, and a small 
value of Vmatch meaus a population of galaxies is only present 
in a small fraction of the volume in at least one of the sur- 
veys. In other words, the larger this value for a CMASS 
galaxy, the more likely this galaxy is a progenitor of a typi- 
cal LRG galaxy, and vice-versa. 

Fig. [S] shows a mapping of the average value of this 
weight onto the two CMASS targeting parameter spaces: 
ag — rvsr — i plot, and a d± vs the cmodel magnitude 
in the i-band. We show the results using the FSPS models 
in the solid contours and the results using Mil in the line 
contours, which are qualitatively similar. The colour-colour 
plot shows a clear trend for the average value of Vmatch to 
increase to redder g — r colours, as expected if LRGs were 
exclusively made of metal rich and old stars. Interestingly, 
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Figure 7. Rest-frame, k+e corrected colour-magnitude relations for CMASS galaxies (filled contours) and LRGs (overploted black 
contours), as a function of Mio.55 shown on the left panel, and as a function of MrO.55 on the right hand side. CMASS galaxies show 
a broader range in their rest-frame ~ ^iO.SSt as well as fainter reach and median in both magnitudes. The right-hand side plot 

shows a clear trend of rest-frame colour with MrO.55, with redder colour going with lower luminosity. This is trend is a result of target 
selection, particularly the magnitude cut - we show the expected evolution of the colour-magnitude relation of an object at the faint end 
of the survey (cmodel= 19.9 at z = 0.45) and an observed colour of r — i = 0.8, between z = 0.25 and z = 0.7. Any object to the faint 
side of the red lines would fail the magnitude cuts of the CMASS algorithm. 



2.0 h J ' I 2.01 ♦ 




g — r cmodel 1 



Figure 8. Average Vmatch weight as a function of colours and i— band magnitude, shown for the two main targeting parameter space 
diagrams in CMASS. A darker colour corresponds to a lower value of Vniatch 1 a-nd the brighter colours to the regions in parameter space 
that have the largest likelihood of being progenitors o f the LRG sample. T he red solid lines show targeting cuts. The orange line on the 
plot on the left shows the morpho logy cut derived in [Masters et al" I 1 I2OIII) . and the dashed green line shows the blue cut of the cut-II 
selection in lEisenstein et al.l 1I2OOII ). 



we also see that some blue regions of the colour-colour plot 
display an increase of the average value of the Vmatch weight. 
This relation is a result of the small but significant amounts 
of young to intermediate-aged stars detected in LRG spectra 
at BOSS redshifts (corresponding roughly to stars aged be- 
tween 1 and 3 Gyr in SDSS-I/II galaxies). The orange line 
in the left-hand plot o f Fig. |S] shows the g — i = 2.35 cut 
of [Masters et al.l (|201ll ). which was motivated by the mor- 
phological analysis of a small subsample of CMASS galaxies 
with Hubble Space Telescope (HST) imaging. They suggest 
selecting galaxies with g — i > 2.35 produces a cleaner sample 



of early-type galaxies (90%) that are more traditionally as- 
sociated with typical LRGs. Additionally, we predict that at 
least a fraction of the galaxies that sit in the blue end of that 
colour-colour plot are also LRG progenitors, temporarily vis- 
iting the blue cloud due to small amounts of star formation. 
Assuming they retain their morphology (it is hard to imag- 
ine a scenario where they would not), our analysis makes 
quantitative predictions on the fraction of star-forming el- 
lipticals that should be found on that part of the diagram, 
given the morphological mixing of the LRG sample (not cur- 
rently known, to our knowledge). This result can be turned 
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Figure 9. k+e corrected absolute magnitudes for CMASS galax- 
ies (black), LRGs (green) and the subset of CMASS galaxies that 
is seen in less than 5% of the LRG volume according to our model 
(purple). These lie almost exclusively at the faint end, demon- 
strating how important the apparent magnitude cut is in the 
sample matching between the two surveys. Solid lines for results 
using the FSPS models and dashed lines for results using Mil. 



into a test of SPS models, as different sets of models will pre- 
dict a different number density at those colours. We leave 
this exploration for future work. 

The right-hand side panel of Fig. [8] shows an uninter- 
rupted trend to lower VInatch towards fainter magnitudes. 
Interestingly, the slope of the Knatch contours are almost 
parallel to the sliding cut in d± with i— band magnitude. 
This cut was designed to follow a line of constant stellar 
mass (Maraston et al. in prep) suggesting that the Knatch 
has a clear dependence on stellar mass, as it should. 

A complementary way to examine the V^match weights is 
to isolate the CMASS galaxies with a small Vmatch weight 
- these are the CMASS galaxies that are less likely to be 
the progenitors of a typical LRG. Fig. |9] shows the K-l-e- 
corrected absolute magnitude distribution of those CMASS 
galaxies with a Vmatch < 0.05, i.e. that are observed in less 
than 5% of the volumes of the surveys. We clearly see these 
galaxies are well confined to the faint end of the CMASS 
population. The difference between the two models is a con- 
sequence of the steeper luminosity evolution given by K-l-e 
corrections of the Mil models - CMASS galaxies are typi- 
cally brighter at LRG redshifts (when compared to a flatter 
luminosity evolution), and are seen through more of its vol- 
ume. 

Fig. [TO] presents the distribution of the absolute rest- 
frame ro.55 — io.55 colour for the same populations as in 
Fig. [5] The bias towards losing intrinsically redder galax- 
ies is explained by the fact that the CMASS sample is itself 
biased towards redder galaxies in Mr-o.55 — Mio.55 at the faint 
end (see Section [3.31 and Fig.[7| due to the i-band selection. 

We show the fraction of CMASS galaxies that are ob- 
served in less than 5% of the LRG volume as a func- 
tion of redshift, absolute magnitude, g — r and rest-frame 
AfrO.55 — Mio.55 colours in Fig. 1111 Once again these fig- 
ures demonstrate that magnitude is the dominant reason 
why these galaxies are not well matched between samples. 



Figure 10. Distribution of k-|-c corrected, absolute ro.55 - *0.55 
colours for CMASS galaxies (black), LRGs (green) and the subset 
of CMASS galaxies (purple) that is seen in less than 5% of the 
LRG volume according to our model. Solid line for results using 
the FSPS models and dashed line for results using Mil. 



but rest-frame colour also plays a part - see the upturn in 
the fraction of lost objects for bright Mro.55 compared to 
the fraction of lost objects for bright Mio.55. These are the 
galaxies with redder Mro.55 — Mio.55 rest-frame colours. 



5 MEASURING POPULATION EVOLUTION 

In order to compute merger and luminosity growth rates, we 
first define the samples of CMASS galaxies and LRGs to be 
investigated (Section l5.1|l . Having selected matched samples, 
we then study the evolution of a number of quantities. In 
Section [5.21 we consider luminosity functions and in Section 
I5.3l the rates of change in number density, luminosity density, 
and typical luminosity per object. 

5.1 Sample selection 

In each survey we take the brightest objects until we reach a 
given K+e corrected absolute magnitude, and we compute a 
Vmatch - weighted comoving number density n and a Fmatch 
- weighted luminosity density £. We consider the following 
options to define the limiting magnitude in each sample, 

Mmin,CMASS aud Mmin,LRG'- 

• A flat cut in k-f E corrected absolute magnitude across 
the two surveys: in this case Mmin,CMASS ~ Mmin,LRG- In 
general, ulrg / ticmass and £lrg 7^ Icmass- 

• A cut in K+e corrected absolute magnitude such that 
both samples have the same comoving number density. In 
this case ulrg = ncMASS by construction, but in general 

Mmin,CMASS / Mmin,LRG , aud (.LRG / (-CMASS- 

• A cut in K+e corrected absolute magnitude such that 
both samples have the same comoving luminosity density. 
In this case Ilrg = £cmass by construction, but in general 
Mmin,CMASS 7^ Mmin,LRG, and TiLRG / TicMASS- This cau 
be advantageous in clustering analyses that are luminosity 
weighted (see Section [S]). 
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Figure 11. The fraction of CMASS galaxies that is seen in less 
than 5% of the LRG volume as a function of redshift (first panel), 
absolute magnitude (second panel - solid line for Mro.55 and 
dashed line for Mio.55, observed g — r colour (third panel), and 
k+E correct rest-frame colour M^o.ss — JV/io.ss (bottom panel). 
Red lines for results using the FSPS models and blue for Mil. 



To avoid confusion we will refer to the number and lumi- 
nosity densities computed using a flat cut in absolute mag- 
nitude as n' and I' . 



5.2 The luminosity function 

With full knowledge of the completeness of the sample, we 
can compute luminosity functions and study their evolution. 
The completeness, in terms of the sample one intended to 
select, is primarily affected by the following well-understood 
effects: 

(i) targeting completeness - not all objects that pass the 
targeting cuts are targeted due to bright star masks, fibre 
collisions or other tiling issues; 

(ii) redshift failure - not all objects with a spectrum suc- 
cessfully yield a redshift; 



(iii) star/galaxy separation - galaxies that fail the star- 
galaxy separation in spite of being genuine galaxy targets. 

We use the targeting co mpleteness and redsh ift fail- 
ure corrections as described in IPercival et al] l|2007t ) for the 
LRGs and in iRoss et al.l (|2012l ') for CMASS galaxies: both 
samples have very high spectroscopic completeness (> 97%). 
The fraction of galaxies lost to the star/galaxy separation 
can be estimated from commissioning data, where star- 
galaxy cuts are less restrictive or not included at all. This 
fraction is estimated to be 1% for CMASS galaxies (Pad- 
manabhan et al. in prep), 1% for cut-I I LRGs and << 1% 
for cut-I LRGs feisenstein et al.l I2OOI] '). This could result 
in a systematic underestimate of the number density of 
CMASS galaxies compared to LRGs, which would a t most 
be ~ 1%. In an independent analyses, [Masters et al.l (|201ll ) 
found 3% ± 2% of CMASS targets in the COSMOS field 
that failed the star-galaxy cuts, in spite of being obviously 
galaxies when captured in high-quality HST imaging. This 
measurement agrees well with the numbers cited above. 

Fig. [T5] presents luminosity functions weighted by 
Vmatch (right), and by the standard V/Vmax weights (left). 
For reference, in both panels we show in the dashed lines the 
luminosity function without any completeness correction - in 
this case it is simply the number count of galaxies per magni- 
tude bin, divided by the volume of each survey. We compute 
the luminosity function in M^o.ss (top) and Mro.55 (bottom) 
absolute magnitudes. Recall that the Vmatch scheme weighs 
each sample such that populations are matched in terms of 
volume, but that it does not yield true volume densities (see 
Section [4. Compared to V^/Vmax weights this downweights 
faint galaxies in both samples, such that the overall lumi- 
nosity functions are matched. In case of zero merger evolu- 
tion or contamination (and in the case of perfect modelling) , 
our Knatch weights fully account for changes in the stellar 
evolution and the two luminosity functions should therefore 
match. Differences can be interpreted in a number of ways: 

(i) growth (i.e., merging); 

(ii) contamination: galaxies in CMASS that have identi- 
cal colour and magnitudes to LRG progenitors but evolve to 
be something else at low redshift; 

(iii) resolution issues: close pairs of galaxies failing to be 
resolved in CMASS due to instrumental and atmospheric 
limitations; 

(iv) inadequacies in the modelling - in this case, mostly 
in the slope of the k-l-E corrections. 

It is clear that the luminosity functions of CMASS 
galaxies and LRGs are better matched in Mio.55 than in 
Mro.55- There is a larger uncertainty in the slope of the 
k-|-E corrections in the r— band, as that traces a region of 
the spectrum sensitive to small amounts of star formation 
at Zc — 0.55. Small mismatches in the amount of star for- 
mation at those redshifts between our composite model and 
the true star-formation rate of CMASS galaxies may not be 
enough to down-weight them using our method, but reveal 
themselves in a detailed comparison such as the one we at- 
tempt here. We therefore argue that the i— band luminosity 
is more reliable for the purposes of our analysis, as it is a 
better tracer of overall luminosity, or stellar mass, of the 
galaxy. 
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Figure 12. Vniatch a-nd Vmax weighted luminosity functions in the k+e corrected ro.55— and jo.55 bands (obtained using the FSPS 
composite model), for the CMASS and LRG samples. The dashed lines show the un-weighted luminosity functions. The Vmax weights 
work by mostly up-weighting the fainter galaxies, as can be seen in the two left panels. This typically breaks down for faint galaxies. 
The Vmatch weight, in turn, up- and down-weights galaxies according to their relative presence on the other survey - this can be seen 
in how effectively we down-weight faint galaxies in both surveys to get a luminosity function that is well matched - particularly in the 
*0.55~band. Poisson errors are negligible (~ 1%) except for the brightest or faintest half magnitudes (1 — 10%). See text for further 
discussion. 



Differences in the shape of the luminosity function can 
help identify the reasons for the differences between the two 
samples. We present a more quantitative analysis in the next 
Section, where we construct three estimators to quantify 
differences in the amplitude and shape of the luminosity 
function, but first we look at the effect of using a different 
k-|-E correction model. 

Fig. 1 131 shows the same luminosity functions as Fig. 1121 
but using the absolute k-l-E corrected absolute magnitudes 
obtained using the Mil models. The differences are sub- 
stantial, especially for the LRGs. Note the differences are 
already apparent in the uncorrected (dashed) curves, show- 
ing the reason lies with the computation of the absolute 
magnitudes themselves, and not with the weighting scheme. 
These results are consistent with the steeper k+E correction 
and the magnitude distributions shown in Figs.[B]and[51 The 
effect is primarily due to the k-l-E corrected absolute magni- 
tudes of the LRGs at Zc = 0.55 - they are « 0.3 magnitudes 
fainter than predicted with the flatter FSPS k-l-E correction. 
The differences are larger for the LRG magnitudes simply 
because of our choice of Zc, which minimises the effect of the 
modelling for CMASS galaxies (see Section 13. 2p . 



There is an overall improvement in the matching of all 
Vmatch luminosity functions across the two surveys when us- 
ing only red CMASS galaxies (with g — i > 2.35) for both 
models. This improvement is small, of only a few per cent, 
and is explained by the fact that the Vmatch weights are 
lower for the bluer galaxies, and so they are already being 
down-weighted when using the full sample. 

Contrasting the two weighting schemes we see that the 
standard V/Vmax weights up-weight galaxies at the faint 
end. Bright galaxies are visible in most of the survey and 
therefore incur a small correction. This shifts the break of 
the luminosity function to fainter magnitudes when com- 
pared to the uncorrected curve, but the falling in number 
density after that must not be trusted completely - V/Vmax 
weights get increasingly dominated by poisson error towards 
faint magnitudes (see Section 14. ip . This is visibly the op- 
posite than what happens using the Vmatch weights in the 
opposite panels. 

5.3 Rates of change 

In order to understand the differences seen in Figs. ll2l and ll3l 
we define three estimators to quantify changes as a function 
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Figure 13. Same as Fig. 1121 but using the absolute magnitudes computed with Mil models. 



of magnitude. For a pair of samples matched on luminosity 
density, we define a merger rate as 



riLRG 



1 



ncMASs J At ' 



(28) 



where At is the time, in Gyr, between the mean redshift 
of the two samples (defined such that At > 0) . Similarly, 
for a pair of samples matched by number density we define 
a luminosity growth as 



re 



Ilrg 



1 

At' 



(29) 



These two rates would be exactly a merger rate and a 
luminosity growth in the absence of complications such as 

(i) resolution issues: close pairs of galaxies failing to be 
resolved within instrumental and atmospheric limitations; 

(ii) contamination: galaxies in CMASS not following our 
composite stellar evolution model and evolving into a differ- 
ent region of colour and magnitude space than that of the 
LRGs at low-redshift; 

(iii) loss of light to the intra-cluster medium (ICM) when 
a merging event occurs; and 

(iv) a systematic offset in the computation of the absolute 
magnitudes as a result of the modelling. 

We investigate (i) in Section [5.41 (ii) is an intrinsic lim- 



itation of any methodology without a full understanding of 
the evolution of all galaxy types, (iii) can potentially be 
investigated by using small-scale clustering and a halo oc- 
cupation distribution type of approach, in order to estimate 
the fraction of satellite merging and a fraction of light lost to 
the ICM. We do not perform such an analysis in the present 
paper, but we will show in Section [7] how, when taken to- 
gether, the results we show in this and in the next Section 
(large scale clustering) present a picture that points strongly 
towards a small amount of population growth. To deal with 
(iv), we also define a galaxy growth rate by using our sam- 
ples matched by a fixed k-l-E corrected absolute magnitude 
(see Section [STTjl as 



^'cMASs/^'cMASS 



1 

At 



(30) 



Tg would match the merger rate even in the presence 
of contaminants (assuming the luminosity function of the 
contaminants was the same as the luminosity function of 
the CMASS galaxies). More generally, it can be interpreted 
as a rate of change of luminosity per single object across 
the two surveys. Whereas rjv and are dominated by the 
relative amplitude of the luminosity function between the 
two redshifts, Vg tells us about differences in the shape. 
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Figure 14. The merger rate, per Gyr, computed as per equation 
1281 as a function of the magnitude of the faintest LRG in the 
sample. The black lines shows rjv x 100 for the full sample, and 
the red line for galaxies with g — i > 2.35. The results obtained 
from using Mil models (dashed lines) show the same slope with 
magnitude as the results using FSPS models (solid lines), but are 
a factor of two to three lower. Poisson errors shown. 
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Figure 15. The luminosity growth, per Gyr, computed as per 
equation 1291 as a function of the magnitude of the faintest LRG 
in the sample. The black lines shows r£ X 100 for the full sample, 
and the red line for galaxies with g— i > 2.35. The results obtained 
from using Mil models (dashed lines) show the same slope with 
magnitude as the results using FSPS models (solid lines), but are 
a factor of two to three lower. Poisson errors shown. 



5. 3. 1 Results 

We compute rjv and as a function of Mm. 55 (the mag- 
nitude of the faintest LRG in the sample, which was used 
to compute the matched samples - see Section [5.1[) . which 
are shown in Figs. [14] and 1151 Our most inclusive samples 
(i.e., where Af^o.ss = —22) include ~ 95% of the LRGs and 
~ 40% of CMASS galaxies, and have large stellar masses 
with logio M/Mq > 11.2 (Maraston et al. 2012, in prep). 
tm is negative for all magnitudes, although it tends 



to zero towards brighter magnitudes. This implies that, for 
the same integrated luminosity density, there are more LRG 
galaxies per comoving volume than there are CMASS galax- 
ies. Le., CMASS galaxies appear to be brighter than LRGs 
in the jo.55 band. This is expected from our analysis of the 
luminosity functions of Figs. ll2l and [T51 We emphasise that 
if this brightening was due simply to the stellar evolution, 
and in the absence of other complications, then our model 
and yniatch weights would account for it. 

ri naturally tells a similar tale - for the same comov- 
ing number density, LRGs are hold less luminosity than 
CMASS galaxies. Removing galaxies with observed colour 
g — i < 2.35 reduces this number by < 1% at the faintest 
magnitudes, but a 5% discrepancy remains, even for the red- 
dest galaxies in the CMASS sample. As is obvious from the 
luminosity functions in Figs. ll2l and ll3l these rates are heav- 
ily dependent on the slope of the k+E corrections. Results 
using the Mil models are identical in shape, but are lower 
by a factor of two to three. Le. - the uncertainty in the 
modelling of the k+E corrections can potentially overwhelm 
these statistics. We return to this at the end of this section. 
One point of interest is how the Vmatch Mio.55 CMASS lu- 
minosity function seems offset from that of the LRGs by an 
almost constant factor as a function of magnitude for both 
FSPS and Mil - this is likely a result of a k+E correction 
slope that is too steep. 

To help understand the observed evolution, we exam- 
ine the rate of change in weighted luminosity per object, or 
rg as given by equation pOf) . which we show on the left- 
most panel of Fig. 1161 Recall that, for this statistic, we se- 
lect galaxy samples based on a fixed k-f E corrected absolute 
magnitude. Using either SPS model, rg is between —1% (at 
the bright end) and 2% (at the faint end). A steeper evo- 
lution seen with Mil is now clear, and it indicates that 
the typical luminosity per galaxy increases between the two 
surveys, especially at the faint end. A similar trend is seen 
using FSPS models, but it is less significant. Processes like 
merging would act to change the shape of the luminosity 
function, according to the fraction and magnitude of the 
merging galaxies. However, that is not what is observed in 
the Afio.55 luminosity function with either set of models. In 
other words, the fact that we observe a small value of rg is 
support for a slowly evolving weighted luminosity per galaxy 
between the two surveys. Note that the sign is positive - i.e. 
I "■LRG^'^LFiQ < 1, or in other words there is on aver- 

"CMASS/'^CMASS 

age more luminosity per galaxy in the LRG sample. This is 
now consistent with a small amount of luminosity growth 
through merging. 

For comparison, we also show rg computed using un- 
weighted number and luminosity densities, or using Vmax 
-corrected densities. In the unweighted case, we see a much 
steeper trend in inferred merger rate with luminosity. This 
trend is dominated by incompleteness issues within the LRG 
sample, which becomes serious at around Afio.55 ~ —23, as 
can be seen in the dashed lines of Figs. ll2l and [THl A 1// Vmax 
weight results in a lower inferred merger rate down to lower 
magnitudes (Mio.55 = —22.5), but shows a steep trend of 
increasing rg with decreasing luminosity beyond that. It is 
difficult to assess whether this effect is due to V/Vmax being 
insufficient to fully correct for completeness or whether it is 
due to a steeper merging rate at those luminosities (which in 
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Figure 16. The change in average Ught per luminosity, per Gyr, computed as per equation 1301 as a function of faintest galaxy in the 
sample. The three panels show the different weighting schemes used when computing number and luminosity densities: Vmatch on the 
left, unweighted in the centre, and Vmax on the right. The black lines shows rg X 100 for the full sample, and the red line for galaxies 
with g — i > 2.35. When weighted by Knatch i shows evidence for a slowly evolving population using bot h stellar population models 
(Mil in dashed lines, FSPS in the solid lines; we also show the purely passive model of lMaraston et al.ll2009l in the dot-dashed red line - 
see Section 15.5.11 for details). The trend in the middle panel is dominated by incompleteness issues in the LRG sample, which are severe 
for MiQ.55 > —23 (see Figs. [121 and I13( l. V/Vmax weights (right) result in a low Vg down to lower magnitudes than V/Vmax , but it rises 
a steeply with decreasing luminosity beyond that. This could be a result of an inadequate completeness correction, or increased merging 
rate at these luminosities. In any case, this comparison demonstrates clearly that the way in which the Vmatch weights balance the two 
samples at low luminosities results in a well matched sample in terms of comoving densities and average luminosity per galaxy - as is 
our goal. Poisson errors are shown for one of the sets of models only for clarity - they are identical for the other set. See text for further 
discussion. 



turn are down- weighted using the VInatch approach). In any 
case, this comparison demonstrates quite clearly that the 
way in which the Knatch weights balance the two samples at 
low luminosities results in a well matched sample in terms 
of comoving densities and average luminosity per galaxy. 

To summarise: we have a complicated scenario: rjv and 
re only reflect a true merger rate or luminosity growth in the 
absence of contamination or unresolved pairs, and a true 
contamination/unresolved pairs fraction in the absence of 
merging. These two quantities are also sensitive to a change 
in the slope of k+E corrections as they rely on matching 
samples by luminosity and number density. They show a 
significant excess of luminosity in CMASS, with respect to 
what we should expect from LRGs. Vg, measuring the change 
in the average luminosity per object, is less sensitive both to 
the slope of the k-l-E correction and to contaminants (pro- 
vided they have a similar luminosity than the galaxies of 
interest). This quantity shows a modest evolution between 
the two surveys (< 2%) for both stellar population synthesis 
(SPS) models. 

We take this investigation further by seeing whether 
unresolved pairs in CMASS could explain the excess of lu- 
minosity implied hy vn and ri alone. 

5.4 Unresolved pairs 

We investigate this issue by looking at pairs of LRGs 
with another object (photometrically classified as a galaxy), 
within a « 2" separation - the angular size subtended by 
same the physical distance at z = 0.3 that corresponds to 
1.2" a.t z — 0.55. In other words, we find all LRGs with a 
close companion such that they would be likely unresolved 



due to seeing (taken to be typically 1.2") at CMASS red- 
shifts. In order to increase our statistics, and to allow us 
to investigate this issue to fainter magnitudes, we perform 
this analysis in the LOZ sample. The LOZ targeting is very 
similar to that of the LRGs in terms of colour, but targets 
fainter galaxies (see Section [2|. We use the full photometric 
sample as this sample is very pure, with stellar contami- 
nation at less that 2% (Padmanabhan et al. in prep). We 
apply the cuts desc ribed in Section [2] on DR8 photometry 
lAihara et al.l [2OIII ). resulting in approximately 1 million 
targets. Of these, only « 15,000 (30,000), or roughly 1.3% 
(2.4%) have a pair between 1.2" and 2" (2.4"). Approxi- 
mately half of these close neighbours are photometrically 
classified as galaxies and half are photometrically classified 
as stars. We also note that the foreground volume of a 2" 
arcsec disc at z = 0.3 is roughly one third of the foreground 
volume of a 1.2" arcsec disc at z = 0.55. So assuming a 
constant number density of foreground objects, we should 
multiply our estimate of the multiple fraction due to chance 
alignment by a factor of 3. We have no way to estimate how 
many of the the close pairs are chance alignments and how 
many are physically associated pairs. If we assume that all 
pairs are chance alignments we reach an estimate on the 
number of unresolved targets at CMASS redshifts of ~ 2%. 
We show the distance profile of these pairs in Fig. 1171 

We can use our estimate of unresolved multiples to cal- 
culate the additional flux brought into the CMASS sample 
from potentially unresolved neighbours. Fig. [18] shows the 
excess in the r— band as a function of r,nod from pairs clas- 
sified as galaxies (black dots, with the median in green). 
We over plot the unnormalised distributions of Vmodei be- 
fore (blue) and after (red) adding the flux of the close pair. 
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Figure 17. Number of LOZ galaxies with a photometric pair as 
a function of its distance. The blue line shows the distances dis- 
tribution for photometric pairs classified as stars, and the red line 
shows the distances for photometric pairs classified as galaxies. 
The black like is the sum of the two. The vertical dashed lines are 
representative of the seeing discs at ^ = 0.3(1.2") and the angular 
size of the 1.2" seeing disc a.t z = 0.6 redshifted to z = 0.3(2"). 
The clear drop off in the number of pairs at distances smaller 
than roughly 1" is due to the fact that we cannot resolve pairs 
closer than the seeing disc. 




Figure 18. The excess magnitude in r— band introduced by po- 
tentially unresolved pairs from neighbours photometrically clas- 
sified as galaxies, as a function of r— band magnitude (black dots, 
with the median in green) . The histograms show the distribution 
of rmod before (blue) and after (red) adding the flux of the close 
neighbour (histograms are not normalised to the y-axis, but share 
a common normalisation). 



It is clear that the effect can be quite dramatic (approxi- 
mately 0.2 magnitudes) for objects fainter than Vmod ~ 18 
magnitudes. It is worth pointing out that in spite having 
of smaller statistics, if we repeat the above analysis around 
LRG targets (as opposed to LOZ galaxies) we get perfectly 
consistent results. 

Integrated over all galaxies with close neighbours, the 
total flux brought into the sample by neighbours photomet- 
rically classified as galaxies is roughly 7.5% of the r— band 
flux of the LRGs. However, this is only happening to roughly 
2% of the CMASS sample according to our more generous 
estimate, and is therefore too small to explain the observed 
excess in luminosity in CMASS, compared to what is ex- 
pected from LRGs if we attribute this excess of luminosity 
to unresolved targets. Reversing the question, to explain the 
excess in luminosity that we observe in the case of FSPS 
models at the faintest end (a 4% excess in luminosity inte- 
grated over the sample; see Fig. [15}, we require that over 
50% of the CMASS galaxies are in fact unresolved targets 
due to chance alignments (this number would have to in- 
crease by approximately a factor of two to explain the ex- 
cess in luminosity inferred using the Mil models). This is 
25 times larger than the fraction estimated by our analyses 
of close pairs in LOZ, suggesting that the slope of k-f-E cor- 
rections or contamination, rather than unresolved targets, 
is the mostly the source for the trends seen in Figs. 1141 and 
1151 and explains why we see only a small evolution in rg in 

[Masters et al.l (|2011l ) identifled a significant number of 
unresolved targets in CMASS by looking at HST COSMOS 
data of a small sub-sample of CMASS galaxies. They show 
that ~ 21 ± 4% of CMASS galaxies are in fact unresolved 
pairs, of which approximately half are estimated to be a re- 
sult of chance alignment, and half physically connected pairs 
(i.e. satellites in the same dark matter halo). Note that we 
are not interested in any unresolved CMASS objects that 
are also unresolved at low redshift, as any such close neigh- 
bours will have their flux accounted for in our estimates of 
the LRGs luminosity, and therefore would not contribute 
towards the discrepancy shown in Figs. [14] and 1151 With- 
out information on the radial distribution of the unresolved 
CMASS targets analysed by Masters et al., a direct compar- 
ison is not particularly insightful. 

Our analysis shows clearly that unresolved targets can- 
not account for the excess in luminosity observed in CMASS 
galaxies. This remains true even allowing for the fraction of 
unresolved pairs measured by Masters et al. Once again we 
emphasise that the slope of k-|-E corrections, rather than 
unresolved targets, is most likely the reason for the trends 
seen in Fie:s. ll4l and ll5l and explains why we see only a small 
evolution in Vg in Fig. 1161 

5.5 Beyond the composite model 

In Section 13.11 we introduced the composite model as the 
best estimate of the overall average colour and magnitude 
evolution of the full LRG sample. In this Section we take the 
opportunity to briefly consider two additional approaches to 
modelling the colours and magnitudes of LRGs. Firstly, we 
consider a purely passive stellar model and secondly we con- 
sider using all individual LRG models instead of averaging 
them into a single composite stellar prescription. 



© 0000 RAS, MNRAS 000, 000-000 



18 Tojeiro et al. 



5.5.1 A passive model 

We ta ke the fully passive stellar model of iMaraston et al] 
l|2009l ) and assume it applies to every LRG and to every red 
CMASS galaxy (with g - i < 2.35). CMASS galaxies bluer 
than this cut show signs of star formation (Thomas et al. 
in prep) and a significantly distinct morpholog ical mix wit h 
a large fraction of late- type galaxies ( Masters et al.|[201l[ ): 
a purely passive model is simply n ot a correct desc r iption 
of the bluer CMASS galaxies. The IMaraston et al.l (|2009l ') 
passive model is based on a single burst of star formation 
at 2; ~ 5, with solar metallicity and an additional compo- 
nent of 3% (by mass) of metal-poor old stars. As this model 
was fitted to 2dF S PSS LRG and Quasar (2SLAQ) data 
l|Cannon et al.ll2006t ). with a significantly redder selection 
function, it does not fit the CMASS data as well as our 
composite model (see Fig. [3|. 

We proceed in exactly the same way as for the compos- 
ite model, and we show the resulting evolution of rg with 
Mio.55 as the dot-dashed red line in Fig. 1161 The new results 
agree particular well with the results using the FSPS com- 
posite model, as expected given Fig. 2] the k-l-E correction 
for a purely passive model follows closely that of the FSPS 
composite model k+E corrections. This shows that the dif- 
ferences in the modelling stem mostly from differences in the 
assumed star-formation histories, but note that these in turn 
are driven by the different stellar evolution tracks assumed 
in each set of stellar population synthesis models. Once we 
weight by Vmatch all models give consistent results. The dif- 
ferences are more significant without the Knatch weight, and 
the inferred merger rates are then larger with a purely pas- 
sive model. 



5.5.2 Using 124 stellar evolution models 

Whereas the composite model captures the average stellar 
evolution of the LRGs, there is a significant amount of scat- 
ter around this average especially in the case of the Mil 
models. We attempt to use the full range of individual fits 
to LRGs of different colour, redshift and luminosity as fol- 
lows. To each LRG we assign the correct model according 
to its r — i colour, luminosity and redshift. To each CMASS 
galaxy we assign the closest model in terms of colour: we 
take the predicted g — r and r — i colours of the 124 models 
at the redshift of each CMASS galaxy, and we assign to that 
galaxy the model that sits closest. Note that as photometric 
scatter is larger than the typical distance between different 
models this is an intrinsically noisy process. 

As explained in Section[37l]each stellar evolution model 
has a different scope in redshift. We therefore must change 
the definition of our weights to allow for this fact: we use 
only the volume probed by each stellar evolution model to 
define Vi^nc, V^^f , Vcmass and V^^^^^ in Eqs. ^ and 
(|27p . This is equivalent to splitting into 124 pairs of surveys, 
each with a different predicted evolution and redshift range, 
and considering the result on the combination. We show the 
resulting evolution of Tg in Fig. 1191 

When compared to the values obtained with the com- 
posite model we see a larger difference in r-g in the case of 
Mil models; this is expected given the larger scatter around 
this composite model. Given how noisy the process of as- 
signing a model to each CMASS galaxy can be, it is hard 
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Figure 19. The change in average light per luminosity, per Gyr, 
computed as per equation 1301 as a function of faintest galaxy in 
the sample, using 124 models of stellar evolution as described in 
Section 15.5.21 The black lines shows Tg X 100 for the full sam- 
ple, and the red line for galaxies with g — i > 2.35. Mil results 
are shown in dashed lines, and FSPS results in solid lines. These 
results should be compared to the left-most panel of Fig. 1161 ob- 
tained using the composite model. The results are similar in the 
FSPS case, but steeper for the Mil models. The differences are 
not significant given the errors. 



to interpret this difference, which is not in any case signif- 
icant. Whereas it is obviously desirable to include as much 
information as possible on the stellar evolution of the LRGs, 
photometric scatter makes this method unreliable. One po- 
tential improvement over the composite model is to con- 
struct two or more average models, for sufficiently distinct 
areas of colour-colour or colour-magnitude space. We leave 
such explorations for future work. For the moment we em- 
phasise that the composite model we describe in Section [3. II 
is the most suitable choice for the current analysis. 



6 LARGE-SCALE CLUSTERING 

The evolution of the large scale clustering has the potential 
to help us interpret the results of the previous sections by 
constraining the merging history of the sample. The evolu- 
tion of the la rge-scale li near bias has a well-defined evolu- 
tion, given bv lFrvl (Il996l ). for pure passive evolution (i.e., no 
mergers). This gives us the opportunity to check whether 
our weighted samples are consistent with a small amount of 
merging and contam ination, a s suggested by F ig llSI 

We will foUow iToieiro fc Percivail (|2010l ') and weigh 
galaxies by their luminosity. The advantage is that any 
merging happening amongst CMASS galaxies will not con- 
tribute towards a deviation from the Fry et al. model, pro- 
vided that no significant loss of light happens to the ICM. 
Note that even if this loss is significant, it is still preferen- 
tial to match samples by luminosity density and weight by 
galaxy luminosity - see Section [T] 
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6.1 Measuring and modelling the correlation 
function 

The two-point correlation function, ^(r) measures the excess 
probability, dP(r) of finding a pair of galaxies at a given 
distance r, compared to a purely random distribution: 
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and follow the same procedure. To ensure a stable inversion 
of the covariance matrix, and to increase our signal-to-noise 
in each bin, we re-bin the correlation functions to 11 bins in 
comoving distance, logarithmically spaced between 30 and 
200 Mpc/h. This results in a total of 22 measurements to be 
fitted by two parameters, totalling 20 degrees of freedom. 



dP{r) = n[l + ^{r)]dV 



(31) 



In practice, we c o unt p airs of galaxies in bins of r and /i, 
and use the Xjmdj^&i^ziyiiy (! 19931) estimator as 

where DD, DR and RR are normalised galaxy-galaxy, 
galaxy-random and random-random pair counts in bins of 
r and fj, respectively is the cosine of the angle between 
a galaxy pair and the line of sight). We use a random cat- 
alogue with the same angular mask as the data catalogue, 
and with an n{z) matched to that of the data. To avoid con- 
tributions from shot noise from the random pair counts, we 
use random catalogues with 10 times the number density of 
the data. 

Setting ^ = in equation (|32|l gives us the monopole 
of the correlation function, as defined in iHamiltonl (|l992l ) 
- this is the excess of finding a pair of galaxies at given 
distance r averaged over pairs observed at all angles with 
respect to the line of sight. The quadropole, or ^ = 2 con- 
tains the next order of information, by effectively comparing 
the power along and across lines of sight, and ^2 are both 
affected by redshift-space distortions and enhanced cluster- 
ing along the line -of-s ight, which we model. Even though the 
passive model of iFrvl l |l996 l constrains only the spherically 
averaged power, or ^o, we fit our data to models of and 
^2, as this improves our signal. 

We m odel the isotropi c , ^-av eraged correlation function 
^(r) as in ISamushia et al] (|2012l ). A non-trivial survey ge- 
ometry imprints a non-uniform distribution of pairs in /i on 
the data as not all galaxy-pair configurations are allowed 
by the window function. We correct for this effect as in 
ISamushia et al.l (|2012t ). by weighting each galaxy pair such 
that the weighted distribution of pairs in /i corresponds to 
that expected in the absence of a window function. 



6.2 Fitting the correlation function 

To increase our resolution with redshift, we split each of 
the CMASS and LRG slices into two, giving a total of 4 
luminosity-matched slices centred at z =0.3, 0.4, 0.5 and 
0.6. For each of the slices we compute ^o('') and (,2{r) accord- 
ing to Equation 1321 and we use a simple 2-dimensional 
minimisation in order to find the best fitting scale-invariant 
amplitudes. 

We estimate the errors and their covariance by using 
mock simulations. We use mock catalogues constructed us- 
ing the Large Suite of Dark Matter Simulations (LasDamas, 
McBride et al. in prep) in order to construct 80 indepen- 
dent realisations of and ^2 for the first two redshift slices 
(we sub-sample each mock in order to reproduce the n{z) in 
each redshift slice). For the last two redshift slices, we us e 
600 perturbation-theory halo mocks of lManera et al.l l|2012l ). 



6.3 Large-scale bias evolution 

The evolution of the amplitude of the monopole can be seen 
in the filled circles of Fig. 1201 with error bars derived from 
the fits to all of the mocks, using the covariance matrices 
described in the previous section. 

To ch eck whether our results are consistent with the 

iFrvl (|l996l ) evolution, we model the red shift evolution o f the 
amplitude of the monopole, Aq{z), as (iHamiltonll 19921 ) : 



Ao{z) 



{z)+'^f{z)b{z) + \f{z)]al{z) (33) 



with as{z) = (Jii{Q)D{z)/D(Q) and 



b{z) = [h{zo) - 1] 



D{zo) 



+ 1 



(34) 



/ is the logarithmic derivative of the linear growth factor 
D{z) with expansion, / = dlog D(2:)/dlog a. For simplic- 
ity we assume a LCDM cosmolo gy and f(z) = with 
7(2) = 0.557 - 0.02z (|Polarski fc Gannouii|[20oi '). 

We take zq = 0.3 and use a simple minimisation to 
fit the model of Equation [33] to our four data points. We 
perform this analysis eight times: 

(i) by weighting the galaxies by their luminosity and 

Vinatch , 

(ii) by not weighting the galaxies; 

(iii) by weighting the galaxies by their luminosity and 
VInatch and applying a. g — i > 2.35 cut; and 

(iv) by not weighting the galaxies and applying a. g — i > 
2.35 cut 

with each of the two stellar population models. The val- 
ues of our fits can be seen in Fig. 1201 For both the FSPS 
and Mil results, the best-fit comes from when the data is 
weighted by Knatch and by luminosity, and when we use only 
galaxies redder than gr — i = 2.35. This result is a good in- 
dication that our weights are performing as expected, and 
that there is less evolution towards the red end of the galaxy 
population, as expected. Weighting by Vmatch and luminosity 
makes a larger difference to the best-fitting than cutting 
the galaxies in colour - this is because the Vmatch weights 
effectively down-weight blue galaxies most effectively. The 
weights increase the overall amplitude of ^0 simply because 
we are up- weighti ng the most luminous objec ts and these are 
more biased (e.g. IZehavi et al.ll2005bl . 120111 ). FSPS models 
give a formally better fit than Mil models, but note that 
in the case of weighted red galaxies, both models give ac- 
ceptable fits to the passive model. This result is a welcome 
confirmation of our interpretation of Fig. 1161 In summary, 
our weights and sample matching yield a sample of galaxies 
that is consistent with dynamical passive evolution. More- 
over, this is robust to the set of stellar population models 
used to create the samples and compute the weights. 
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Figure 20. The amplitude of the luminosity and Vmatch -weighted large-scale amplitude of go computed as a function of redshift (solid 
black dots). The two lowest redshift points lie exclusively within the LRG sample, and the two highest redshift points exclusively in the 
CMASS sample. The open dots show the amplitude of the un-weighted large-scale power, fitted to the same scales. The red dots (open 
and fillod) show the amplitude of go when selecting only CMASS galaxies with g — i > 2.35. The lines show the best fit passive model of 
IPrvl (1996) (Equation I33I I. obtained by assuming LCDM and fitting for b(zQ) - the solid line is a fit to the filled (weighted) points, and 
the dashed line is a fit to the open (not weighted) points. We give the minimum values of for each case. Left: FSPS models; right; 
Mil models. 



7 SUMMARY, DISCUSSION AND 
CONCLUSIONS 

In this paper we present a joint analysis of SDSS-I/II LRGs 
and BOSS CMASS galaxies, with the aim of identifying and 
characterising a coeval population of galaxies spanning a 
redshift range between 0.23 and 0.7. We are motivated from 
the desire to select a population of galaxies that is evolving 
passively in a dynamical sense (i.e. no mergers) as closely 
as possible, as the large-scale bias evolution of such a pop- 
ulation can be understood analytically and yield significant 
gains in cosmological analysis of large-scale structure. 

We focused on the progenitors of LRGs, as massive red 
galaxies are prime candidates for such a population. As the 
targeting selection in CMASS is significantly wider both in 
colour and magnitude (Section [2]), we developed a set of 
weights that optimally matches galaxies in terms of their 
stellar evolution (Section I4.1|l . To do so, we relied on the 
fossil record of LRGs from which we extracted a stellar evo- 
lution model (Fig. E]). We then used this model to iden- 
tify the most likely progenitors of LRGs amongst CMASS 
galaxies (Section l4.2|l . Finally we developed a number of es- 
timators to attempt to characterise and quantify population 
evolution between the two surveys (Sections [5] and |6]). 

We find that CMASS galaxies are an extension of LRGs 
by being both intrinsically fainter and having a wider range 
in observed and in rest-frame colour. Fig. [TT]shows the frac- 
tion of CMASS galaxies that are not expected to evolve 
into LRGs as a function of redshift, absolute magnitude, 
observed g — r colour and rest-frame Mro.55 — Mio.55 k-l-E 
corrected colour. We find a steep dependences in absolute 
magnitude and rest-frame colours, confirming that CMASS 
galaxies are broader in terms of intrinsic properties. 

Our analysis of weighted number and luminosity den- 
sities using the rjv and ri estimators (Figs. [T^] and I14p 



points towards a scenario where the CMASS sample is typ- 
ically brighter expected from the LRG progenitors. To in- 
vestigate this issue further we considered the potential con- 
tamination from unresolved targets - i.e. CMASS targets 
that are targeted as a single object, but are in fact unre- 
solved pairs of stars of galaxies (Section 15. 4|l . By examining 
close pairs of LOZ galaxies (such that they would be un- 
resolved at CMASS redshifts) we estimate that the extra 
luminosity from these pairs would be too small (< 1% in 
the r"'^^— band) to explain the excess in luminosity we see 
in CMASS, even in our most generous scenario. The most 
likely reason for the remaining differences is uncertainty in 
the slope of the k-l-E corrections (which affect the absolute 
magnitudes mostly of the LRGs), or contaminants - objects 
in CMASS that share a region of colour-magnitude space 
with LRG progenitors but that evolve into something other 
than present-day massive red galaxies. Even though both 
stellar population synthesis models used in this study give 
an identical trend of this luminosity excess with magnitude, 
we find it to be 2-3 times larger with the Mil models. To 
help identify the reason for this luminosity excess, we look 
at the rate of change in luminosity per object, r^, and at the 
evolution of the large-scale clustering. 

The estimator rg was designed to be intrinsically less 
sensitive to contaminants in CMASS (provided they have a 
similar luminosity distribution to the galaxies of interest, in 
which case rg can be interpreted as a merger rate) . It is also 
less sensitive the slope of the k-|-E corrections (mostly via 
the sample selection, which is different to the sample selec- 
tion needed for rjv and re - see Section [5TT}. We find that rg 
tests differences in shape of the luminosity function, whilst 
rjv and re are also sensitive to the relative amplitudes. We 
only small evolution in rg. Moreover, we find that this evo- 
lution (between —1% at the bright end and 2% at the faint 
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end - see left-hand panel of Fig. I16|l is much less sensitive 
to the stellar population synthesis modelling. We compared 
the optimally Knatch -weighted Vg with the results obtained 
using a standard F/Vmax weight (Fig. [T6)l . In that case we 
would infer a merger rate of up to 13% at the faintest end, 
showing that the weighting scheme we introduce is having 
the effect we intended: effectively matching the two sam- 
ples terms of comoving densities and average luminosity per 
galaxy. This result is evidence for a small evolution in the 
properties of the galaxies, once they are properly weighted. 

We should comment on a possible systematic error in rg 
due to the completeness estimation in each survey. In Sec- 
tion [5]2] we noted that the fraction of CM ASS galaxies that 
fail the star-galaxy separation (but which are nonetheless 
genuine galaxies) is larger than the fraction of LRGs that 
fail similar cuts at low redshift by at most 1%. This is a 
small number, but comparable to the accuracy with which 
we aim to constrain the dynamical evolution of these galax- 
ies. We should therefore conservatively add 1% uncertainty 
to the rates we show in Section 15.3.11 For our most inclu- 
sive sample (Mio.55 < —22), rg is therefore constrained to 
be 1.6 ± 1.5% with Mil models, and 0.4 ± 1.4% with FSPS 
models. Our most inclusive samples include « 95% of the 
LRGs and ~ 40% of CMASS galaxies, and have large stellar 
masses with logn, M/Mq > 11.2 (Maraston et al. 2012, in 
prep). 

To place further constraints on the evolution of our 
weighted samples, we investigate the evolution of the am- 
plitude of the large-scale clustering (Fig. l2Up . We find that 
the best fit to a dynamically passive evolution model (i.e. 
strictly no mergers, and assuming LCDM) happens when 
we weight the galaxies by Knatch and luminosity, and we use 
only red galaxies {g — i> 2.35). This is further evidence that 
the weights are working as they should, and we find formally 
good fits to the passive model using both FSPS and Mil so- 
lutions, providing increased weight to our interpretation of 
Figs. [151 [14] and [H 

7.1 Comparison with previous work 

Independent efforts to constrain the assembly and evolu- 
tion of massive galaxies range generally focus on the cosmic 
evolution of the luminosity function or on their clustering. 
Direct comparisons are difficult due to the different sam- 
ples used, which vary in terms of number density (or typical 
minimum halo mass), magnitude or colour cuts, but this is 
non etheless a u s eful e xercise . 

ICool et al.l l|2008l ) measured the evolution in the lumi- 
nosity function of LRGs up to z « 0.9 by observing a small 
sample of nearly 300 galaxies targeted using similar cuts to 
SDSS-II LRGs. They find that at the massive end, galaxies 
must not have grown by more than 50% between z = 0.9 
and z = 0.1 , corres ponding to less than 8% growth per Gyr. 
IWake et all (|2006^ analysed SDSS-I/II and 2SLAQ LRGs, 
covering a redshift range between z — 0.17 and z = 0.6 and 
found that no growth was necessary to explain the evolution 
of the luminosity function once a passive stellar evolution 
model was applied, but that non-passive growth with up to 
25% of the galaxi es merging between t he two redshifts could 
not be ruled out. iBrown et al.l (|2007l l studied the luminos- 
ity evolution of nearly 40,000 galaxies since z ~ 1 in the 
Bootes field. At the massive red end, they find that at least 



80% of the stellar mass of these galaxies had to be place 
by z = 0.9, resulting i n a mild growth of l ess th an 3% per 
Gyr to the present-day. iToieiro fc Percival (I2OI0I ) measured 
the stellar growth of LRGs between 0.2 < z < 0.45 and 
found it to vary between 2-6%, depending on the luminos- 
ity with brighter galaxies showing a smaller evolution. For 
the same luminosity range in this paper we find a smaller 
merger rate, in spite of the larger redshift range. This dis- 
crepancy is explained by a combination of different model 
for the stellar evolution of the galaxies and by an estima- 
tor that is less sensitive to overall offsets in magnitudes and 
potential contaminants (see equation [30}. Whereas studying 
the evolution of luminosity densities can certainly test de- 
parture from a purely or nearly passive dynamical evolution, 
it only formally tests a net influx of luminosity in or out of 
the samples between the two redshifts. Studies of the evolu- 
tion of the clustering of galaxies are therefore important in 
characterising the mass assembly of these massive galaxies. 

By fitting a halo model to a sample at high-redshift, 
evolving it passively to a lower redshift, and comparing 
it directly to a model fitted to data at that redshift al- 
lows for a direct comparison of predicted (from a passive 
model) and measured populations of centrals and satellites 
for samples matched in number density. Traditionally such 
approaches find too many satellites at low redshift, result- 
ing in too large a number d ensity and an overe stimation of 
powe r at all scales (see e.g. IWhite et a l. 2007; Brown et al.l 
l2008l : IWake et al.ll2008D . 

The usual interpretation is that a fraction of the satel- 
lites and/or centrals must have merged. Such approaches 
have implied merger rates in the order of 2-8% per Gyr since 
a redshift of 1, which is in apparent agreement with the re- 
sults presented here. An interesting disagreement, however, 
is the fact that such studies have found that the evolution 
of the large-scale power departs significantly from the pas- 
sive model of i Frvl (Il996l ): they find a nearly invariant large- 
scale amplitude as a function of redshift, which corresponds 
to effectively underestimating the bias at low redshift, with 
respect to passive evolution. As mentioned previously the 
merging scenario - coming from HOD fits to smaller scales 
- suggests that a fraction of galaxies in the sample at high- 
redshift must have merged, effectively reducing the large- 
scale power via two mechanisms. It primarily reduces the 
number of objects in high-mass halos, therefore decreasing 
the overall bias of the sample. A likely secondary effect is 
that merging within the sample must reduce the number 
density at low redshift for the same population of galaxies - 
as the samples are matched in number density, resulting in 
an enrichment of less-biased galaxies at low redshift. 

Our matching and weighting by luminosity bypases 
both of these problems. Firstly, when weighting galaxies by 
their luminosity, merging events between galaxies will not 
decrease the relative contribution of the halos within which 
they reside to the overall bias of the sample. This is only 
strictly true in the case of no loss of light to the ICM, but 
we argue here that weighting by luminosity will almost al- 
ways be better than any weighting scheme that depends on 
the number of objects - in a merging of two objects the rel- 
ative contribution of a given halo will be reduced by 1/2 
if weighting by number. It follows that, provided that the 
overall loss of light is less than 50% of the combined light 
of the merging system, we have an estimation of the bias 
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evolution that is less sensitive to merging of galaxies within 
the sample, and to which the Fry model is more applicable. 
Note that this is true even in the case of merging within 
the sample. Estimates on how much light may be lost to the 
intr a-cluster medium v aries. Based on a halo model analy- 
sis, IWhite et al.l (|2007l ) estimated an upper bound on the 
loss of light to the ICM of 25 - 40% of the light of the 
accreting satellite (no t of the whol e merging system), de- 
pending on halo mass: Iskibba et al] (j2007) argued that this 
fractio n should be between 5 — 15%. Using simulation-based 
models IConrov et all l|2007l ) show how scenarios that allow 
disrupted satellites to deposit up to 80% of their stars in 
the ICM favour the observed ev olution in the galaxy stel- 
lar mass function since z ^ 1. iPurcell et al.l J2OO3) esti- 
mate that around 20% of the total light in massive halos 
(M > 10^'^ M0) is in a diffuse component. Analyses as the 
ones above show how understanding the mechanisms that 
lead to the formation of the ICL is of clear importance to 
learning how massive galaxies assemble. Nonetheless it re- 
mains clear that weighting by luminosity is in practice ad- 
vantageous if one wishes, as we do, to use a passive model 
for the evolution of the large-scale bias. The effectiveness of 
our weighting scheme is nicely demonstrated with our clus- 
tering analysis, which shows how a passive model is a better 
fit to the weighted data (see Fig. I20|) . 

7.2 Final remarks and future work 

We conclude that our sample is slowly evolving (to less that 
2% by merging, when samples are appropriately matched 
and weighted), to the extent of what is testable by current 
data and models. We demonstrated the efficiency of a Knatch 
and luminosity weighting in constructing a sample that is 
as close to being dynamically passive as possible, whilst at 
the same time not needing to cut the sample in colour or 
redshift; therefore optimising our signal. This aspect is par- 
ticularly important given broad nature of CMASS galaxies 
with respect to LRGs. 

In terms of future work, we will assess the cosmological 
gains of using our slowly evolving sample when measuring 
growth rates and redshifts-space distortions using the large- 
scale amplitude of the correlation function. We further in- 
tend to extend this analysis to the LOZ and main galaxy 
samples jStrauss et al.l [20021 ) . which will provide not only 
better statistics for cosmology analysis, but will also allow a 
study of evolutionary paths of other populations of galaxies. 
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